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A new approach is presented that treats NMR spin-spin relaxation as kinetics in spin phase space. 
The approach is applied to free induction decay (FID) in solids containing equivalent nuclear spins 
1/2. The description obtained does not involve adjustable parameters. As an example, the calcula- 
tion is performed for ^®F FID in CaF2, and the results are in good agreement with experiment. 
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We present a new approach to the problem of free in- 
duction decay (FID) in the lattice of equivalent spin 1/2 
nuclei, where nuclear spin-spin interaction is responsible 
for the relaxation. The approach applies the framework 
of the Boltzmann kinetic description beyond the limit 
of instantaneous collisions. The resulting method is for- 
mally comparable to the original treatment of Lowe and 
Norberg [Q, in the sense that it starts from the Ising- 
like Hamiltonian, matches the second and the fourth mo- 
menta, and does not involve adjustments to experimental 
data. Moreover, the T-criterion, which is introduced in 
this work, is also satisfied in the Lowe and Norberg calcu- 
lation. The difference is that our mathematical construc- 
tion is based on physical arguments, which give better 
control over the calculation, and, in principle, allow a 
routine analysis of experiments to be performed to ex- 
tract unknown microscopic parameters. None of the ap- 
proximate methods suggested so far has been generally 
accepted for such a task. 

A typical FID formulation in solids assumes that each 
nucleus is at rest on its site in the crystal lattice, and the 
following conditions are fulfilled: 
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where T is the temperature of the initial equilibrium dis- 
tribution, r2 is the Larmor frequency in the static mag- 
netic field, T2 is the time scale of the nuclear spin-spin 
interaction, and Ti is the spin-lattice relaxation time. 
The spin-lattice relaxation can be neglected, since the 
time range of interest is of the order of T2. 

We consider FID in the Larmor rotating reference 
frame, where the Hamiltonian is H 
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Ik is the spin 1/2 operator of the fcth nucleus, and Akn 
and B^n are the interaction coefficients. The z-axis is 
chosen along the direction of the Larmor precession. 

The FID is measured as a spin response to 7r/2 radio 
frequency pulse. The tt/2 pulse rotates the equilibrium 
spin system so that immediately after the pulse, all spins 
are uniformly polarized in the direction perpendicular to 



the static magnetic field. We choose the a;-axis along the 
magnetization at this moment of time. After the pulse, 
the a;-component of the magnetization decays as a 
result of the spin-spin interactions. The leading term in 
the high temperature expansion gives jl]j~ 
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where 7 is the gyromagnetic ratio. Usually, FID is pre- 
sented as a normalized function G{t) = Mx{t)/Mx{0)- 

Since all spins are equivalent, each of them equally 
contributes to the magnetization. For the purposes of 
the forthcoming consideration, we express G{t) in terms 
of the average magnetic moment m of one spin: 
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We outline our approach by comparing it with the 
derivation of the Boltzmann kinetic equation (BKE) 

The BKE can be considered as a modification of the 
exact equation for the system of noninteracting classical 
particles. In our method, we modify the exactly soluble 
case of the Ising-like Hamiltonian. The BKE is derived by 
counting particles entering and leaving small volumes of 
the phase space. Similarly, we analyze the averaged spin 
behaviour in a given configuration of neighbors. The con- 
figuration of neighbors refers to the phase space domain 
arising in the context of the Ising-like Hamiltonian. 

The derivation of BKE is based on the approxima- 
tion of intantaneous collisions, which is not applicable 
to the nuclear spin-spin interaction in solids. Instead, we 
substantiate the quantitative claim of our approach by 
adopting the criteria presented below. 

The first criterion is that the mathematical construc- 
tion of our theory has to be time reversible. This allows 
a well-defined correspondence to be established between 
the parameters of the theory and the time reversible 
microscopic dynamics. In particular, it enables us to 
match the exactly calculated momenta AI2 = — ^(^L=o 
and M4 = ^j^L=o: which can be obtained by the direct 
trace evaluation of the first terms in the time expansion 
of Eq. (^ . Time reversibility implies abandoning simple 
first-order differential rate equations. Consequently, the 
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minimal description has to be based on time-reversible 
second-order differential rate equations. 

Among the various dynamical correlations that com- 
plicate the analysis, the most important are the two-spin 
correlations. If an approximate calculation can accu- 
rately take the two-spin correlations into account, then 
it is reasonable to expect that the contribution from the 
higher order correlations is more random, i.e better av- 
erageable. The criterion, which at least partially guar- 
antees that the two-spin correlations are respected, is as 
follows: The FID shape obtained by the effective calcu- 
lation is the function of the microscopic coefficients A^n 
and Bkn- We require that if in this function all Akn and 
Bkn, except for the coefficients describing the interac- 
tion between any given pair of spins, equal zero, then the 
function reproduces the exact FID shape of the two-spin 
system 
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We call this the "T-criterion." 

The FID can be evaluated in a closed form in the Ising- 
like case when all Bkn — 0. We use the Lowe and Norberg 
interpretation Q of this evaluation to initially motivate 
the formalism of our description. 

In the Ising-like Hamiltonian, the operator of the local 
field, which affects the kth spin, is 
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Since the z-component of each spin is the constant of 
motion, Eq.(^) implies that each spin is rotated by a 
constant local field created by spin's neighbors. 

If all operators Im are diagonal in the basis chosen for 
the trace evaluation in Eq.(^), the contributions to the 
trace can be interpreted as a result of spin precessions in 
classical local fields. The possible values of these fields 
are the eigenvalues of the field operator hk ■ 

We introduce index C to refer to the configurations of 
neighbors, and define configuration as a particular set of 
the eigenvalues of operators 1^ in Eq.(||). The descrip- 
tion can be restricted to a finite number of neighbors Nq. 
Consequently, there are infinitely many spins surrounded 
by a given configuration C of Nq neighbors. In the fol- 
lowing we use notation to indicate averaging over 
all configurations. 

For configuration C, we define mcx{t) as the average 
magnetic moment of nuclei that are surrounded by neigh- 
bors with the specified set of instantaneous spin projec- 
tions on the z-axis. The initial uniform polarization im- 
plies that mcx{0) = m^{0) = 

The evolution of mcx is governed by the equation 
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Each configuration of neighbors uniquely specifies the 
combination of signs in Eq.(^). 

As a result, mcx{t) = rncx(0)cos{ujct) , and Eq.(|4|) with 
mx{t) = {mcx{t))c leads to 
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The right-hand side (RHS) of Eq.Q is the Fourier 
transform of the discrete distribution of ojc, which can 
be considered as a sum of random quantities ±>lfe„/(27i). 
Therefore, according to the central limit theorem, the 
distribution approaches Gaussian when 
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(5) where M2 = ^^^En ^L- In this case, the RHS of Eq.(|) 



can be approximated as exp{—M2 t /2). 

Apart from the Ising-like case, the evolutions ofmcx in 
different configurations of neighbors are mutually depen- 
dent. We employ the following generalization of Eq.(0) 
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where Pcx+ and Pcx- arc the probabilities that the spin 
projections on the a;-axis are positive and negative, re- 
spectively; Wc+ is the effective rate of the transition from 
a state where the spin is oriented positively along the x- 
axis to the state where the spin is oriented negatively 
along the same axis; Wc- is the rate of the reverse tran- 
sition; mo = 7?i/2 is the maximum magnetic moment of 
one spin. The probabilities obey the relationships: 



Pcx+ + Pcx- - 1, 

'mo{Pcx+ - Pcx-) = mcx 



where ojc is the precession rate in configuration C: 
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The rates in Eq. (|11|) characterize the net effect of two 
factors: the direct influence of neighbors and the flux of 
spin polarization to and from conflguration C. When the 
spin system is weakly polarized, each of these factors can 
lead to a shght asymmetry between Wc+ and Wc-- 

The first factor can be understood after we rewrite the 
two-spin Hamiltonian extracted from Eq.(^ as 

Wfcn = Bknikxinx ~ '^Akn{Ik+In+ + Ik- In-) 

+ ^{Akn + 2Bkn){Ik+In- + Ik-In+), (14) 

where I-^- = ly + il^^ , and = ly — il^ ■ The second and 
the third terms in Eq.(p^ lead to the double-fiip of two 
parallel spins and the flip-flop of two antiparallel spins, 
respectively. If the fcth spin is parallel to the average 
spin polarization, it is more probable that the nth spin 
will be parallel to the fcth spin. This effectively increases 
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the influence of the double-flip term on the kth spin and 
reduces the influence of the flip-flop term. If the fcth spin 
is antiparallel to the average polarization, the opposite 
effect would occur. When rUx > 0, the double-flip corre- 
lation increases W^_^_ and reduces in Eg. (pi]). The 
flip-flop correlation leads to the opposite result. 

The other factor mentioned above is frequently referred 
to as motional narrowing. In terms of our description, 
this means that the faster the time variations of the real 
local fields, the greater the difference between the cho- 
sen configuration of neighbors and those configurations 
that actually drove the spin to its current surroundings. 
One can conclude that spins coming to configuration C 
from other configurations tend to be slightly polarized in 
the direction of the average magnetization. This reduces 
W^_(_ and increases W^_, provided rrix > 0. 

We define the average rate as 
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The high temperature condition guarantees that the dif- 
ference between and W^_ is small. The basic as- 
sumption of our analysis is that the leading term of this 
difference can be expressed as 
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where a is a configuration independent parameter, which 
is determined by the average of all factors considered. 
The status of assumption (^^ is somewhat similar to 
the self-consistent relaxation time approximation, which 
is frequently adopted to solve BKE. 

The substitution of Eqs.(|l|,|l|,|l|,|l6l) into Eq.(|l] 
yields 
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Given the initial uniformly polarized state, the solution 
of Eq.(|l|) is 

mcx{t) = mcx{0) cos{Wct) 

- a I mx{t - t')Wcsin(Wct')dt' . (18) 
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Smcemx{t) = {mcx{i))c^ the averaging of Eq. (p^ over 
all configurations, together with Eq.(^), results in the 
integral equation 



G{t) = g{t) + a fcit - t')^4^dt' 
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where g{t) — {cos{Wct))i,. 

We choose g{t) and a such that the second and the 
fourth momenta obtained from Eq.(|l^) match the exact 
calculation. Matching the second moment gives 
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where Mag = -|lL=o - 
moment, combined with Eq.(E 
pressed as 
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where M4, = ^L^„ = (I¥4)c 



(20) 

Matching the fourth 
allows a to be ex- 
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Eqs.( pO| , plD relate the shape and the scale of the ac- 
ceptable distributions of Wc ■ Namely, the shape specifies 
the ratio of M^g/M^g, which enters Eq.( [2l|) a nd deter- 
mines the value of a. With a known a, Eq.(|2C|) gives the 
value of A/2g, which defines the scale of the distribution. 

An important property of Eqs.(20 2^) is that they 
guarantee the fulfillment of the T-criterion for any shape 
of the distribution of Wc with finite moments. The exact 
solution (H) for the system of two spins 1/2 always gives 
the ratio M4/M2 = 1, which formally leads to an infinite 
value of a in Eq. ( ^l]). The divergence does not appear in 
the solution of Eq.(|ig|), because, according to Eq. (po|) , it 
is offset by the small value of M2g. When the two-spin 
problem is considered as a limit of the many-spin problem 
with Akm Bkn 0, except for A12 and B12, the solution 
of Eq.(|l^) converges to the RHS of Eq.(||), independently 
of the shape of the distribution of Wc ■ 

The relevance of the T-criterion is supported by the 
fact that the observable part of G{t) obtained from 
Eqs.( p^ , p0|j2l] ) is weakly sensitive to the variations of the 
input shape of the Wc distribution. Thus relatively crude 
assumptions about this shape still should lead to good 
accuracy in the result. 

Since the fulfillment of the T-criterion is guaranteed, 
and the important many-spin correlations are taken into 
account by the choice of parameter a, we assume that 
the distribution of Wc is produced by uncorrelated con- 
tributions from spin neighbors. As a result, this distri- 
bution has a tendency to be Gaussian. However, in the 
presence of a few strongly interacting neighbors, it can 
have a peaked structure originating from the Ising-like 
part of Hamiltonian (H) . This part alone would produce 
the discrete distribution of the static rates wc given by 
Eq.(||). The discrete structure, presumably, propagates 
to the distribution of Wc , but it has to be "washed out" 
by the polarization flux between different configurations. 
The proper scale this "lifetime" effect is given by 
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where f{t) = Tr{e^ ^ * /ife e"^ " * hk}. 

We choose the input shape of the Wc distribution to be 
the convolution of the distribution of luc with a Gaussian 
that has the second moment M2/. Based on this shape. 
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g{t) becomes a renormalized product of a Gaussian and 
the RHS of Eq.(g): 

g{t) = exp[-^M2f{vtf] Hcos {^^A^nVt) , (23) 



where rj ■ 



is the renormahzation factor 



(1 + Q)(M2+Af2/) 

required by Eq.(pO|). If the large number of neighbors 
criterion (p^ is satisfied, g(t) becomes Gaussian. 

Eqs. ( |19| , |2l| , p^ , p3| ) form a closed system, which allows 
G{t) to be calculated. 

We apply the method to CaF2, where FID is observed 
on ^^F nuclei ( 7 — 25166.2 rad s~^ Oe~^), which form 
a simple cubic lattice having near-neighbor separation of 
d = 2.72325A. The spin-spin interaction is assumed to 
be magnetic dipolar 1^. 

The calculation is performed with the static magnetic 
field oriented along the [100], [110] and [111 ] cr ystal di- 
rections. Evaluating g{t) according to Eq.(|23|), we in- 
dividually include 50 cosines with the interaction con- 
stants making the largest contribution to M2g, and ap- 
proximate the product of other cosines by the Gaussian 
exponent, which adds the remaining contribution to M2g. 
We also approximate M2 / by the second moment of one 
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is sufficient given the alternating signs of A^n- This gives 
M2/ — 1^2. Variations of M2/, even by a factor of 2, 
result in a negligible difference in the observable part of 
FID. In Eq.@, we take the exact ratios M^/M^ from 
Ref. §, and obtain the values of a: (0.54, 0.42, 0.42) for 
the [100], [110] and [111] directions, respectively. 

When g{t) and a are known, Eq.(^9|) is easily soluble 
numerically. The solutions are plotted in Fig. [T], together 
with the experimental data of Engelsberg and Lowe [^ . 

The microscopic coefficients originating from the mag- 
netic dipolar interaction have the relation Bkn — ^\Akn- 
As a result, the double-flip correlation, which is dis- 
cussed after Eq.(^4|), is sufhcicntly strong to outweigh the 
flip-flop correlation and the motional narrowing. Conse- 
quently, a > 0, which leads to the oscillating FID shape. 

The value of a can also be negative, e.g. when Akn 
and Bkn have the same sign for each pair of spins. If 
g{t) is Gaussian, and a changes from to —1, the shape 
of G{t) obtained from Eq. (p^) evolves from Gaussian to 
nearly exponential. The minimum possible value of a 
is —1, which leads to G{t) = 1 for any g{t). This limit 
arises as the interaction approaches the Heisenberg form 
with all Akn = 0. 

In summary, we presented a kinetic FID theory, which 
reproduced the exact results of the two-spin problem with 
Hamiltonian and many-spin problem with Ising-like 
and Heisenberg Hamiltonians. Based on the physical ar- 
guments, the theory extended the above exactly solu- 
ble cases in the space of all possible interaction coeffi- 
cients Akn and Bkn- Given also the fact that the first 
four derivatives of G{t) at i = obtained from Eq.(p^) 
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FIG. 1. FID in CaF2 with static magnetic field along [100], 
[110] and [111] crystal directions. Solid lines are the solutions 
of Eq.(p^. Dashed lines are experimental data from Ref. j^. 

coincide with the exact calculation, we conclude that in 
most cases, the solution of Eq.(|l9|) approximates the ob- 
servable part of FID with good accuracy. The accuracy 
is mainly limited by the adequacy of assumption (|l^) . 
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